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Abstract 

We present a 1.91457-approximation algorithm for the prize-collecting travelling salesman 
problem. This is obtained by combining a randomized variant of a rounding algorithm of 
Bienstock et al. [2] and a primal-dual algorithm of Goemans and Williamson [5] . 

1 Introduction 

In the prize-collecting travelling salesman problem (PC-TSP), we are given a vertex set V (with 
\V\ = n), a metric c on V x V (i.e. c satisfies (i) = Cji > for all i,j 6 V and (ii) triangle 
inequality: Cij + Cjk > c%k for all i,j,k G V), a special vertex r £ V (the depot), penalties 
7r : V — > M+, and the goal is to find a cycle T with r £ V(T) such that 

c(T)+tt(V\V(T)) 

is minimized, where c(T) = S(i j)eT c *i> 7r (5') = Sies 71 "*) an d ^(^) denotes the vertices spanned 
byT. 

The first constant approximation algorithm for PC-TSP was given by Bienstock et al. [2]. It is 
based on rouding the optimum solution to a natural LP relaxation for the problem, and provides a 
performance guarantee of 2.5. Goemans and Williamson [5] have designed a primal-dual algorithm 
based on the same LP relaxation, and this gives a 2-approximation algorithm for the problem. 
In 1998, Goemans [4] has shown that a simple improvement of the algorithm of Bienstock et 
al. gives a guarantee of 2.055- • • = - e _ 2 /3 • Recently, Archer et al. |T) are the first to break the 
barrier of 2 and provide an improvement of the primal-dual algorithm of Goemans and Williamson; 
their performance guarantee is 1.990283. In this note, we show that by combining the rounding 
algorithm of Bienstock et al. and the primal-dual algorithm of Goemans and Williamson, we can 
obtain a guarantee of 1.91456 • • • = 1 _a^_ 1/3 ■ The analysis uses the technique in [4] together with 

an improved analysis of the primal-dual algorithm as observed in [3] and used in Archer et al. pQ . 
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2 Combining Approximation Algorithms 



We start by briefly reviewing the rounding result of Bienstock et al. [2j. Consider a classical LP 
relaxation of PC-TSP: 

Min y^c e x e + y~]ir(v)(l - y v ) 

e£E v 

subject to: 

x{5{v)) = 2y v veV\{r} 
(LP) x{5(S)) >2y v S C V, r $ S, v € S 

< x e < 1 e G E 

0<y v <l v EV 

Vr = 1, 

where E denotes the edge set of the complete graph on E. For conciseness, we use c(x) + 7r(l — y) 
to denote the objective function of this LP. Let x* , y* be an optimum solution of this LP relaxation, 
and let LP = c(x*) + tt(1 — y*) denote its value. Bienstock et al. [2] show the following (based on 
the analysis of Christofides' algorithm due to Wolsey [8] and Shmoys and Williamson [7]). 

Proposition 1 (Bienstock et al.). Let < 7 < 1 and let S('j) = {v : y* > 7}. Let T 7 denote the 
cycle on S^) output by Christofides' algorithm when given S (7) as vertex set. Then: 

c(T 7 ) < ^c(x*). 

The 2.5-approximation algorithm can then be derived by setting 7 = | since we get c(T 3/ / 5 ) < 
|c(x*) and ir(V \S(3/5)) < |7r(l — y*). In [3], we have shown that one can get a better performance 
guarantee by taking the best cycle output over all possible values of 7; notice that this leads to at 
most n — 1 different cycles. 

The primal-dual algorithm in [5] constructs a cycle T and a dual solution to the linear pro- 
gramming relaxation above such that their values are within a factor 2 of each other, showing a 
performance guarantee of 2 since the value of any dual solution is a lower bound on LP. Chudak, 
Roughgarden and Williamson [3] (see their Theorem 2.1) observe that the analysis of [5] actually 
shows a stronger guarantee on the penalty side of the objective function, namely that the cycle T 
returned satisfies: 

c(T) + (2- -^—] n(V \ V(T)) < ( 2 - LP. (1) 



n — 1 J \ n — 1 

This increased factor on the penalty side is exploited in Archer et al. [I], and this motivated the 
result in this note. Suppose now that we apply the primal-dual algorithm to an instance in which 
we replace the penalties n(-) by ir'(-) given by 

Av) = ^-^ ( V Mv). (2) 

2 — l/(n — 1) 

Thus, ([1]) implies that the cycle T returned satisfies: 

c(T) + n(V \ V(T)) < (2 - -L-) LP', (3) 
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where LP' denotes the LP value for the penalties vr'(-). As the optimum solution x*, y* of LP (with 
penalties vr(-)) is feasible for the linear programming relaxation with penalties ir'(-), we derive that 
the cycle T pd output satisfies: 



Summarizing: 

Proposition 2. The primal-dual algorithm applied to an instance with penalties tt'(-) given by ^ 
outputs a cycle T pd such that 



We claim that the best of the algorithms given in Propositions Q] and [2] gives a better than 2 
approximation guarantee for PC-TSP. 

Theorem 3. Let 



where a = - — 2 — rp? < 1.91457. 

l-fe-V3 

As mentioned earlier, the minimum in the theorem involves only n different algorithms as we 
need only to consider values 7 equal to some y*. 

Proof. We construct an appropriate probability distribution over all the algorithms involved such 
that the expected cost of the solution produced is at most a (c(x*) + 7r(l — y*))- 

First, assume that we select 7 randomly (according to a certain distribution to be specified). 
Then, by Proposition [H we have that 




c(T pd ) + vr(y \ V(T pd )) < 2c(x*) + vr(l - y*). 



H = min(min(c(T 7 ) + tt(V \ V(j))), c(T pd ) + ir(V \ V{T pd ))). 



Then 



H < a (c{x*) + tt(1 - y*)) = aLP, 



E[c{T,)\ < \e i 
2 7 



while the expected penalty we have to pay is 



E[tt(V \ V( 7 ))) = Pr[l > V*(v)Hv). 



v&V 
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Thus, the overall expected cost is: 

"1" 

c{x") + 

Assume now that 7 is chosen uniformly between a = e -1 / 3 = 0.71653 • • • and 1. Then, 



£[c(T 7 ) + 7r(F\F( 7 ))]<~£ 



ix*) + ^2Prh>y*(v)]n(v). (4) 



E 

and 





-/ 


.7. 


J a 



1 1 , ln(a) 1 1 

-—ax — — 



Pr[y > y] 
Therefore, (j3|) becomes: 



I -ax ~ I- a 3(1 -a) 3(1 - e' 1 / 3 ) 



1 < lES 0<y<a. 



E[c(T 7 ) + vr(y \ F( 7 ))] < 2(1 _ 1 e _ 1/3) c(x*) + x.g-i/a ^ 1 " ^" ( 5 ) 

Suppose we now select, with probability p, the primal-dual algorithm as given in Proposition 
[2] or, with probability 1 — p, the rounding algorithm with 7 chosen randomly acording to 7 ~ 
{7[e _1 / 3 j ij. From ([5]) and Proposition [21 we get that the expected cost E* of the resulting algorithm 
satisfies: 

E* < (2p + (1 - P) 2{1 _\_ 1/3) ) c(xl + (p + (1 - P)TT^) ^ X " 

Choosing p = (1 ~ p) 2{1 _l-i/ 3 y i.e. p = ^^173 > we g et 

< 3p(c(x*) + tt(1 - y*)) = 3pLP. 

Therefore, the best of the algorithms involved outputs a solution of cost at most 3pLP = aLP 
where 

a = 5 — -. - < 1.91457. 

1 - ge^ 1 /^ 

□ 

One can show that the probability distribution given in the proof is optimal for the purpose of 
this proof; this is left as an exercise for the reader. 

Theorem [3] shows that the linear programming relaxation of PC-TSP has an integrality gap 
bounded by 1.91457; in contrast, the result of Archer et al. [1] does not imply a better than 2 
bound on the integrality gap. 

As a final remark, if we replace Christofides' algorithm with an algorithm for the symmetric 
TSP that outputs a solution within a factor (3 of the standard LP relaxation for the TSP then the 
approach described in this note gives a guarantee of 

1 



1 - ie 1 " 2 /^ 



for PC-TSP. 
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